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Abstract. We calculate the generating functions for the number of 
tilings of rectangles of various widths by the right tromino, the L tetro- 
mino, and the T tetromino. This allows us to place lower bounds on the 
entropy of tilings of the plane by each of these. For the T tetromino, we 
also derive a lower bound from the solution of the Ising model in two 
dimensions. 



1 Introduction 



Tilings of the plane are of interest both to statistical physicists and recreational 
mathematicians. While the number of ways to tile a lattice with dominoes can 
be calculated by expressing it as a determinant [1, 2], telling whether a finite 
collection of shapes can tile the plane at all is undecidable [3, 4]. 

In the 1994 edition of his wonderful book Polyominoes [5], Solomon W. 
Golomb states that the problem of determining how many ways a 4 x n rectangle 
can be tiled by right trominoes "appears to be a challenging problem with rea- 
sonable hope of an attainable solution." Here we solve this problem, and several 
more. Specifically, we find the generating functions for tilings of 4 x n and 5xti 
rectangles by the right tromino, and tilings of 4 x n rectangles by the L and T 
tetrominoes. 

The number of ways of tiling a m x n rectangle with any finite collection 
of shapes, where m is fixed, can be found by calculating the nth power of a 
transfer matrix whose rows and columns correspond to the various interface 
shapes a partial tiling can have. As n grows large, the number of tilings N(n) 
grows as A™ where A is the largest eigenvalue of this matrix. In general, the 
number of interface shapes and therefore the size of the transfer matrix grows 
exponentially in m. However, if a small number of interface shapes suffice to 
generate all possible tilings, then we can use a small transfer matrix, and A will 
then be an algebraic number of low degree. In fact, we find that for the problems 
considered here the largest matrix we need is 4 x 4. 

By calculating A for these rectangles, we can place a lower bound on the 
entropy per site of tilings of the plane by these polyominoes. In the case of the T 
tetromino, we improve this bound by using the exact solution of the Ising model 
in two dimensions. 



2 The right tromino 



Consider tilings of a 4 x n rectangle with right trominoes, where n is a multiple 
of 3. The number of interfaces seems potentially large. However, it turns out 
that we only need to think about three, so that a 3 x 3 matrix will suffice. These 
are shown in Figure 1 , with the various kinds of transitions that can take place 
as the tiling grows to the right. 





Fig. 1. The three kinds of interfaces for the right tromino on rectangles of width 4, 
straight, deep jog, and shallow jog, and the transitions between them. 

For these three interfaces, which we call 'straight,' 'deep jog,' and 'shallow 
jog,' we define N(t), N\(t) and A^fi) respectively as the number of tilings when 
there arc n = 3t columns to the left of the dotted line. Obviously Ni(t) and 
N2(t) stay the same if we count tilings of their vertical reflections instead. Our 
goal is to find N(t). To express these, we write them as generating functions, 

G{z)=Y J N{t)z t 
t 

and similarly for G\{z) and G2(z). 

From Figure 1, we see that a straight interface can make a transition to itself 
in four ways, each of which increases t by 1. It has one transition each to a deep 
jog and its reflection, each of which has one transition back to a straight, and 
both of these increase t by 1. Deep and shallow jogs of a given orientation have 
two transitions in each direction; as we have defined them a shallow-to-deep 
transition increments t but a deep-to-shallow transition does not. A shallow jog 
has two transitions to its reflection, incrementing t. Finally, N(0) = 1 since there 
is exactly one way to tile a 4 x rectangle. 

This gives us a 3 x 3 transfer matrix, and the system of linear equations 



G = 1 + 4zG + 1zG x 



Gi = zG + 2zG 2 
G 2 = 2Gi + 2zG 2 



The solution for G is 



G(z) = 



l-6z 



1-Wz + 22z 2 + 4z 3 
The coefficients N(t) of the first few terms of G's Taylor expansion are 

1, 4, 18, 88, 468, 2672, 16072, 100064, 636368, 4097984, 26579488, . . . 
giving the number of ways to tile rectangles of size 4 x 0, 4 x 3, 4 x 6, and so on. 
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Fig. 2. Fault-free rectangles of right trominoes of width 4 and length 3, 6, 9, and 12. 



To gain a better understanding of this series, we concentrate on fault-free 
rectangles, those whose only straight interfaces are those at the left and right 
end. The first few of these are shown in Figure 2. These always begin and end 
with transitions from a straight edge to a deep jog and back, except for t = 1 
where we have a pair of 2 x 3 rectangles. By disallowing transitions back to a 
straight interface except at the end, we find that the linear equations for fault- 
free rectangles are 

G' = 4z + 2zG; 
G[ = z + 2zG' 2 
C*2 — 2>G^ -\- 2zG"2 



This has the solution 

x ~ 2-llz-2z 2 
& W = 2Z l-6z 
The Taylor series of G' tells us that the number N' of fault-free tilings is 

4,2,8,48,288,1728,... 

which we can write as 

C 4 if t = 1 

AT'(i) = < 2 if i = 2 

[ 8 • 6*- 3 if t > 3 

Since any tiling consists of a concatenation of fault-free tilings, we have 

G(z) = 1 + G'(z) + G'(zf + G'(zf + ■■■= 1 _ 1 GI{z) 

which the reader can verify. Note that we have to define G'(0) = N'(0) = in 
order for this formula to work. 

The asymptotic growth of N(t) is the reciprocal of the radius of convergence 
of G's Taylor series. Thus N(t) oc A* where A is the largest positive root of 



This gives 



where tt/2 < 8 < n and 



A 3 - 10A 2 + 22A + 4 = 



A = \ ( 5 + V34 cos V - 



Numerically, we have 



3 /519 
tan^--^ — 



A = 6.54560770847481152029 • • • 



It would be nice to find a simpler expression for A, perhaps using the decompo- 
sition into fault-free rectangles. 

For 5 x 3t rectangles, four interface shapes suffice. We call these straight, little 
jog, big jog, and slope, and their transitions are shown in Figure 3 (reflections 
and reversals are not shown). If their generating functions are G, Gi, G2 and 
G3 respectively, our 4x4 transfer matrix is 

G = 1 + 8zG 2 + 4zG 3 
Gi = 4zG + zGi + 2z 2 G 2 + (2z + 4z 2 )G 3 
G 2 = 2Gi + zG 2 + 4zG 3 
G 3 = 2zG + 4zGi + (2z + Az 2 )G 2 + Az 2 G 3 



3t 



3t 



3t 
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Fig. 3. The four kinds of interfaces for the right tromino on rectangles of width 5, 
straight, big jog, little jog, and slope, and the transitions between them. 



whose solution is 

1 - 2z - 3lz 2 - 40z 3 - 20z 4 



G(z) 



l-2z- 103z 2 + 280z 3 + 380z 4 
The coefficients N(t) of the first few terms of G's Taylor expansion are 



1, 0, 72, 384, 8544, 76800, 1168512, 12785664, 170678784, 2014648320, 25633231872, . 

Note that there are no 5 x 3 tilings. 

The generating functions for fault-free tilings obey 

G' = 8zG' 2 + 4zG' 3 

G\ =4z + zG\ + 2z 2 G' 2 + (2z + 4z 2 )G' 3 

G' 2 = 2G[ + zG' 2 + 4zG' 3 

G' 3 = 2z + 4zG[ + (2z + 4z 2 )G' 2 + 4z 2 G' 3 



and so 



G'(z) = 24z 2 3 + 10z + 15z 2 

1 ' l-2z- 31z 2 - 40z3 - 20z 4 



The first few N'(t), starting with t = 2, are then 

72, 384, 3360, 21504, 163968, 1136640, 8283648, 58791936, 423121920, 
These are shown for t = 2 and t = 3 in Figure 4. 
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Fig. 4. Fault-free rectangles of right trominoes of width 5 and length 6 and 9. 



The asymptotic growth of the number of fault-free rectangles is N'(i) oc A'* 
where A' is the largest positive root of 

A' 4 - 2A' 3 - 31A' 2 - 40A' - 20 = 

which has a rather complicated closed-form expression which we will not repro- 
duce here. Numerically, 

A' = 7.16235536278185348653 • • • 

For all rectangles, including faulty ones, we have N(t) oc A* where A is the largest 
positive root of 

A 4 - 2A 3 - 103A 2 - 280A - 380 = 



This is 



i v 62? y 

* bye. 
where 



Numerically, we have 



i-3yi- 13107 ^4369 + 418,Uyl + 2304 

1 F r— 1 



3 3/ 

209 y3-y 3 -4369 



2 V3yH 



y = 1204327 - 72 V263697405 



A = 12.36366722455963019234 • • • 



3 L tetrominoes 



We can perform a similar analysis for tilings of 4 x n rectangles with L tetro- 
minoes. In fact, here the analysis is easier. As Figure 5 shows, we have just two 
kinds of interfaces, straight and jogged. We define N(t) and N\(i) as the num- 
ber of tilings for a straight or jogged interface respectively when there are 2t 
columns to the left of the dotted line. Then a straight has two transitions to 
itself that increase t by 1, four to itself that increase t by 2, and two to itself 
that increase t by 3. It has two transitions each to a jog and its reflection, one of 
which increases t by 1 and the other by 2. Finally, a jog has two transitions to 
another (reflected) jog that increase t by 1, one transition to itself that increases 
t by 3, and the two inverse transitions back to a straight. 



2t 



2t 
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Fig. 5. The two kinds of interfaces for the L tetromino on rectangles of width 4, straight 
and jogged, and the transitions between them. 



Thus we have a 2 x 2 transfer matrix, and defining generating functions as 
before gives us the linear equations 

G = 1 + (2z + Az 2 + 2z 3 )G + (2z + 2z 2 )G l 
d = (z + z 2 )G+ (2z + z 3 )G 1 

so 

G(z) '- 2Z - Z3 



1 - Az - 2z 2 + z 3 + Az A + Az 5 + 2z 6 



The coefficients N(t) of the first few terms of G's Taylor expansion are 

1, 2, 10, 42, 182, 790, 3432, 14914, 64814, 281680, 1224182, . . . 

giving the number of ways to tile a rectangle of size 4x0,4x2,4x4, etc. 

As before, we can focus our attention on fault-free rectangles. For t > 3, any 
such rectangle consists of transitions from a straight edge to a jog, and transitions 
between jogs in between. For t = 1, 2, 3 we have an additional 2, 4, and 2 fault- 
free rectangles respectively as shown in Figure 5. The linear equations for the 
fault-free generating functions are thus 

G' = 2z + 4z 2 + 2z 3 + 2(z + z 2 )G[ 
G[=z + z 2 + (2z + z 3 )G[ 

This gives 

2z(l + z) 2 (l-z-z 3 ) 
° {Z) = l-2z-z3 

whose Taylor expansion tells us that the number N' of fault-free tilings is 

2, 6, 10, 18, 38, 84, 186, 410, 904, 1994 . . . 

The first few of these are shown in Figure 6. The reader can verify that G(z) = 
1/(1 -G'{z)). 

The asymptotic growth of the number of fault-free rectangles is N'{t) oc A'* 
where A' is the largest root of 

A' 3 - 2A' 2 -1 = 

which is 

w i / 3/43 3 Vrn 3 /43 zVrn] 

3 y V "2" — 2- + Vt + ^J 

or numerically 

A' = 2.20556943040059031170 • • • 

For all rectangles, including faulty ones, we have N(t) oc A* where A is the largest 
root of 

A 6 - 4A 5 - 2A 4 + A 3 + 4A 2 + 4A + 2 = 
This appears not to have a closed form, but it is approximately 



A = 4.34601641141649282849- 
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Fig. 6. Fault-free rectangles of L tetrominoes of width 4 and length 2, 4, 6, 8, and 10. 




Fig. 7. The two kinds of interfaces for the T-tetromino, straight and jagged, and the 
transitions between them. 
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Fig. 8. The T tetromino has two fault-free rectangles of each size. 



4 T tetrominoes 

Of these three polyominoes, finding the number of tilings of a 4 x n rectangle is 
easiest for the T tetromino. In fact, such tilings only exist if n is a multiple of 4. 
Figure 7 shows two kinds of interfaces, straight and jagged, with At columns to 
the left of the dotted line. A straight edge must make a transition either to a jag 
or its reflection, incrementing t. A jag can make a transition back to straight, 
leaving t fixed, or to itself, incrementing t. Thus the generating functions obey 

G = 1 + 2Gi 

d = zG + zd 

This gives 

G(z) 

whose Taylor expansion tells us that 




if t = 
if t > 



Thus we can find a closed form for N(t), unlike in the previous two cases. There 
are exactly two fault-free rectangles for each t, as shown in Figure 8. Thus the 
number of tilings is 2-3* since there are two choices of initial fault-free rectangle 
and three choices for each increment of t, namely either continuing the current 
fault-free rectangle, or ending it and starting one of two new ones. 

It is easy to show that the T tetromino cannot tile any rectangles of width 
5, 6, or 7. Proofs are shown in Figure 9. More generally, Walkup [6] showed that 
a rectangle can be tiled with T tetrominoes if and only if its length and width 
are both multiples of 4. 



l-z 

1 -3z 



5 Tilings of the plane 

Since the plane can be divided into rows or columns, the entropy per site of 
tilings of rectangles of a given width serves as a lower bound for the entropy per 
site of tilings of the plane. Specifically, if the number of ways to tile a m x n 
rectangle is N(m,n), we define the entropy a as 

In N(m,n) 
a = hm 



Fig. 9. Proofs that the T tetromino cannot tile rectangles of width 5, 6, or 7. Attempted 
tilings always lead to one of the un-tileable notches shown on the left, or to an even 
more obviously un-tileable space. 



so that 



N{m, n) oc e 



for large rectangles. For the right tromino, our analysis of 5 x n rectangles gives 



where the factor of 15 comes from the fact that each increment of t adds 15 sites 
to the rectangle. Similarly, for the L tetromino on 4 x n rectangles, we have 



Better lower bounds could be obtained by looking at wider rectangles. Note that 
if there are any tilings we failed to see, this only improves these lower bounds, 
since additional transitions can only increase A. 

For the T tetromino, the lower bound we get from Ax n rectangles, ^ In 3 w 
0.06866, is not as good as the bound | In 2 w 0.08664 which we can derive by 
noting that the 8-cell shape in Figure 10 can tile the plane, and in turn can be 
tiled in two ways by the T. We can get a better bound as follows. 

Suppose that the boundaries between tiles at points {Ax, Ay) are clockwise 
or counterclockwise fylfots as shown in Figure 11. Whenever neighboring fylfots 
have opposite orientations, there is only one way to tile the space between them, 
but if they have the same orientation the space between them can be tiled in 
two ways as in Figure 10. (Note that this is also the source of two of the choices 
in the 4xn tilings above.) If we define the two orientations as +1 and —1, then 
we can sum over all configurations of fylfots {si = ±1}, with each one giving us 
1 choice for pairs of unlike neighbors and 2 choices for like ones. Thus we have 
a lower bound of 



a > — In A = 0.1676508- • • 



cr> -lnA = 0.183657- •• 

8 




Fig. 10. Since this eight-cell shape can tile the plane and in turn can be tiled in two 
ways, the entropy of tilings of the plane with T tetrominoes is at least | In 2. 
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Fig. 11. The space between two fylfots 4 sites apart can be tiled in 1 way if they are 
of opposite type, and 2 ways if they are the same. 



where the product is over all pairs of nearest neighbors, and the sum is over the 
configurations of a lattice of size ^ x j = We can rewrite this as 

N(m,n)= Yl ne^ s ^ +1 ) ln2 

{s» = ±l} ij 

= 2W II e ( 5 ln2 ) s ^ 

{s i= ±l} ij 

(note that there are ^ edges in the fylfot lattice). 

We now note that the latter sum is the partition function of an antifcrro- 
magnetic Ising model with (3 = ^ In 2. We can transform this to a ferromagnetic 
model by negating the Sj on one of the checkerboard sublattices. Using the exact 
solution of the Ising model in two dimensions [7] , we then have 

O > (In 2 + (Thing) 

where 

1 1 f 2vr f 27r 

cising = In 2 + -- — / / dwi dw2 In [cosh 2 2/3 — (sinh 2/3) (cos ui\ + cos u-z)] 

2 (2ny J J 

For /3 = \ In 2, where cosh 2 2/3 = f| and sinh 2(3 = f , we have 



o-ising = 0.8270269567- 



and so 



cr > 0.09501088358- •• 



For upper bounds, we can generalize an argument given in [2] as follows. 
If we construct a tiling of the plane with right trominoes by scanning from 
top to bottom and left to right, at each step the first unoccupied site can be 
given a tromino with only four different orientations. Since we have at most 
n/3 such choices, the entropy is at most |hi4 w 0.462. Similar considerations 
give j In 8 w 0.520 and jln4 w 0.347 for the L and T tetromino respectively. 
Obviously the gap between our upper and lower bounds is still quite large. 

For dominoes or 'dimcrs,' the entropy is known to be G/tt = 0.29156 • • • where 
G is Catalan's constant [1, 2]. It is tempting to think that other polyominoes 
might have exact solutions; one related model which does is the covering of 
the triangular lattice with triangular trimers [8]. However, while the general 
problem of covering an arbitrary graph with dimers can be solved in polynomial 
time, covering it with triangles is NP-complcte [9] , so generalized versions of the 
tromino problem are probably hard. 
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Note added. Meyerowitz [10] has calculated the generating function for the 
T tetromino in rectangles of width 8. 
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